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In order to expand the study on flow instability of supercritical circulating fluidized bed (CFB) boiler, a new nu- 
merical computational model considering the heat storage of the tube wall metal was presented in this paper. The 
lumped parameter method was proposed for wall temperature calculation and the single channel model was 
adopted for the analysis of flow instability. Based on the time-domain method, a new numerical computational 
program suitable for the analysis of flow instability in the water wall of supercritical CFB boiler with annular 
furnace was established. To verify the code, calculation results were respectively compared with data of commer- 
cial software. According to the comparisons, the new code was proved to be reasonable and accurate for practical 
engineering application in analysis of flow instability. Based on the new program, the flow instability of super- 
critical CFB boiler with annular furnace was simulated by time-domain method. When 1.2 times heat load distur- 
bance was applied on the loop, results showed that the inlet flow rate, outlet flow rate and wall temperature fluc- 
tuated with time eventually remained at constant values, suggesting that the hydrodynamic flow was stable. The 


results also showed that in the case of considering the heat storage, the flow in the water wall is easier to return to 


stable state than without considering heat storage. 
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Introduction 


The combustion technology of circulating fluidized 
bed (CFB) has the advantages of low NOx emission, 
wide adaptability, high combustion efficiency and low 
cost of pollution control, which has been widely used in 
many fields, such as electric power, petroleum, chemical 
and waste disposal and so on. Compared with subcritical 
CFB boiler, the operating parameters and mode of super- 
critical CFB boiler are high, and the fluid in the tube is 
likely to operate at the supercritical state of high load, or 
it may work at the two phase state of lower load. In order 


Nomenclature 

A inner cross section area, m? Zz 
dn inner diameter, m Az 
f friction factor 


to ensure the safe and reliable operation of the equipment, 
we must ensure that the film boiling(DNB)doesn't occur, 
and the temperature of tube is in the safe range after the 
drying up (DRYOUT). 

In recent decades, the flow instability has caused wide 
concern on the design and operation of various heat 
transfer equipments, such as nuclear reactors, refrigera- 
tion plants and boilers. It may result in oscillation of wall 
temperature, lead to fatigue damage of tube and heat 
transfer deterioration [1]. Thus making investigations on 
flow instability is very essential. Since 1938, the domes- 
tic and foreign scholars have carried out extensive and 


axial length, m 


space step, m 


Greek symbols 


g gravity, m/s? p density, kg/m? 
h specie Sathaioy: Tike 0 ae en direction with respect to horizontal 
i space index Og dimensional Dirac delta function, m” 
j time index a convective heat transfer coefficient (W/m*-K) 
k pressure drop coefficient Subscripts 
L total length of channel, m in inlet 
M mass flow rate, kg/s out outlet 
n total number of control volume jb local 
P pressure, Pa n inner 
AP pressure drop from inlet to outlet, Pa f fluid 
qi linear power density, W/m w wall 
F: internal surface area per unit length Superscripts 
Cw specific heat(W/m?-K) 0 the initial steady value 
My mass per unit length(kg/m) j the old time value 
t temperature( °C) j+1 the new time value 
At time step, s 


in-depth study on the two-phase flow instability [2-16]. 
However, there is few research on the problem of the 
flow instability considering the heat storage of the tube 
wall metal in a boiler. The heat flux density of the fur- 
nace is influenced by the heat storage of metal, then the 
dynamic characteristics of the furnace will also be influ- 
enced, so it is necessary to study the hydrodynamic flow 
instability of a boiler when the heat storage of the tube 
wall metal is considered. In this paper, the numerical 
computational model considering the heat storage of 
metal and suitable for the analysis of flow instability and 
the calculation of wall temperature is established by the 
lumped parameter method. Based on Fortran, the single 
channel model and the time-domain method are adopted. 
In order to provide guidance for the design and operation 
of boiler, we select the loop that is most likely to occur 
flow instability in the water wall as the research object. 


Mathematic model 


Generally, there are mainly two kinds of approaches 
for solving the governing equations of mass, momentum 
and energy. The first one is frequency-domain method, 
which converts the linearized equations to transfer func- 
tions by Laplace transform [17-23]. However, the fre- 
quency domain method can’t solve nonlinear problems 
very well. The second one is time-domain method, based 
on numerical discretization and integration, which con- 
serves the nonlinear information of the original equations 
[24-28]. The time domain method is employed in the 
present code. The temperature of the tube wall can be 
calculated by using lumped parameter method or distri- 


buted parameter method. In this paper, the lumped para- 
meter method is used to analyze and calculate the wall 
temperature. 


Basic assumptions 


For the purpose of analyzing the flow instability of 
boiler, the following assumptions are proposed: 

(1) One-dimensional model is employed with consi- 
dering compressibility and thermal expansion. 

(2) The temperature and velocity of water are un- 
iformly distributed in cross section, and water flows only 
along the axial direction. 

(3) Only heat transfer in the radial direction is consi- 
dered. 

(4) Effects of kinetic, potential energy and viscous 
dissipation on energy equation are ignored. 


Transient equations for supercritical fluid 


Supercritical water has no phase change which can be 
regarded as single phase. The basic governing equations 
of supercritical water are the same as the equations of 
homogeneous flow model. 

The mass equation is: 

OM JOP aig (1) 
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The momentum equation is: 
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The energy equation is: 
O(Mh) +A O(ph) _ a (3) 
Oz Ot 
The state equation is: 
= f(h,P) (4) 
The metal heat storage equation is: 
Ot 
Hw, (5) 
qi =h = Cy ôt 
The heat transfer equation of working fluid side: 
qn =aF(ty-tr) (6) 


Numerical calculation method 


Discretization 

As shown in Fig1, the channel is divided into several 
control volumes from the inlet to the outlet with equal 
length of Az. The non-staggered grid method is used for 
dividing the control volume of flow channel. The fluid 
physical state parameters, such as mass flow rate, pres- 
sure, density, temperature and enthalpy are located in the 
center of the control volume. Discretization of the go- 
verning equations employs first-order upwind scheme in 
space phase and implicit scheme in time phase. Accord- 
ing to the first-order upwind scheme, physical state pa- 
rameters on the junction are equal to that of the adjacent 
upstream control volume. Integration of the governing 
equations on every control volume yields the following 
equations: 
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Fig. 1 Schematic diagram of control volume 


Mass equation: 


j+1 j 1 j+1 
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Momentum equation: 
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Energy equation: 
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State equation is: 
pi = fhf, P) (10) 
The metal heat storage equation is: 
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The heat m ees of working fluid side: 
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Facilitate the heat transfer equation: 
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Boundary conditions 


For the single channel model, the boundary condi- 
tions are as follows: 

(1)inlet enthalpy hj, is given and be constant; 

(2)inlet pressure Pin is given and be constant; 

(3)pressure drop from inlet to outlet AP is constant. 


Solution method 


A summary of the solution method is shown in Fig. 
2: 

(1) The initial values of P, M, h, p, tw, tp are firstly 
computed. After that, a small perturbation of heat flux 
is imposed on the steady state and steps (2)—(5) are 
repeated to calculate the transient values of the next 
time step level. 

(2) Combine the equation of metal heat storage (11) 


and heat transfer (13), +/*! of the control volume and 


wi 


qi; are solved. 


(3) Assume the inlet mass flow rate M/*'. 


(4) Repeated steps (a) ~ (e) from inlet to ee which 
means i is from 1 to n. 
jti 


(a) Assume the value of o7} . 


(b) According to the mass equation(7) , M/*! is 


i+l 

solved. 

(c) According to the momentum equation(8) , pin 
is solved. 

(d) According to the energy equation(9) , hii is 
solved. 

(e) According to the state equation(10) , an improved 
value of p!ġ! is obtained by using PH, hj! . Com- 


pare the density of the former and latter. Repeat steps 


(b)-(e) until the absolute error of density reaches a de- 
rj+l 


sired convergence. p;/} can be treated as an updated 


guess. 
(5) After Step (4), the outlet pressure P/*' is ob- 


n+l 


tained. 


If the relative difference (2 -Pİ 3)- AP) / AP meets 


in n+l 
the requirement of precision, then the assumption of 
M/*' is correct. Solution turns to next time step level. 


If not, a chord secant method is used to iterate on M/*! 


until the relative difference is satisfied. 
The solution method for steady-state equations is sim- 
ilar to the transient equations which are not mentioned here. 


Validation 


To validate the present code, calculation results ob- 
tained by the present numerical code are compared 
with the simulation results of Dynastab software. Dy- 
nastab software is a kind of commercial software de- 
veloped by Siemens Ltd for simulating the hydrody- 
namic flow, aiming at providing guidance for the de- 
sign and operation of boiler. In this paper, the geome- 
tric structure and thermal parameters are the same with 
those of Dynastab software in practical engineering 
application of a 600 MW supercritical W-shaped boiler. 
Initial calculating parameters are shown in Table.1. In 
order to simulate the flow instability precisely and 
calculate the wall temperature within a smaller error 
range, the factors of inclined angle, heat distribution 
and type of tube are considered, comparison results are 
shown in Fig.3, Fig.4, Fig.5. 

Fig.3 shows that the wave curve of the inlet flow 
versus time derived from the present program is con- 
sistent with the curve derived from Dynastab, it indi- 
cates that the method of this paper has certain reliabil- 
ity for simulating the flow instability of the supercriti- 
cal water. 

According to the Fig.4, we can calculate the relative 
error of the fluid temperature and the wall temperature 
between program and Dynastab is 0.04% and 0.2%. It pro- 
ves that the present program is reasonable to simulate the 
fluid temperature and wall temperature within error limits. 


The initial values are firstly computed. 
Perturbation of heat flux is imposed on the steady state. 
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Combine the equation of metal heat storage 
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Fig.2 Program calculation process 


Table 1 Initial calculating parameters 
Parameter Value 
Inlet pressure 12.566 MPa 
Inlet enthalpy 1256.6 kJ/kg 
Inlet mass flow rate 0.1367 kg/s 
Inlet pressure drop coefficient 0.5 
Outlet pressure drop coefficient 1 
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Fig. 3 The inlet flow versus time 


Fig.5 shows that the wave curve of the inlet flow ver- 
sus time in the case of considering heat storage of the 
tube wall metal is quicker to restore to the initial state 


than the flow fluctuation without considering heat sto- 
rage. It means that the flow in the water wall is more 
stable with heat storage. This phenomenon indicates that 
the heat storage of metal is an influence on the dynamic 
characteristics of the furnace, so it's necessary to consider 
the heat storage of metal when analyze the flow instabil- 
ity of a boiler. 
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Fig.4 The outlet temperature distribution of section 
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Fig.5 Effect of heat storage on the inlet flow 


Analysis of flow instability in supercritical CFB 
boiler with annular furnace 


Annular furnace 


Fig.6 shows the simple structure of the annular fur- 
nace. The furnace of boiler is divided into inner and out 
ring, the outer and inner ring are arranged in parallel. In 
front wall of the outer ring of furnace, three groups of 
cyclone separator are set, the same arrangement in back 
wall. Not only more water cooling walls can be arranged 
in this structure to obtain enough heat exchange areas, 
but also this structure reduces the height of furnace and 
reduces the power consumption of the circulating pump. 

In this paper, we selected the sixth loop of the front 
wall as object of study, the initial calculating parameters 
are expressed in Table.2. The tube in the loop is vertical 
smooth tubes with low mass flow, of which the inner and 
outer diameters are 19mm and 32mm, while the length is 
55.6026m. 
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Fig. 6 Structure of annular furnace 


Fig.7 and Fig.8 show the pipeline structure and heat 
load distribution, the distribution of heat load along the 
height of the furnace is not uniform. According to the 
characteristics of heat flux distribution along the furnace 
height, the water wall is divided into more sections where 
the heat flux experiences a sharp heat flux change and is 
divided into fewer sections where the heat flux expe- 
riences a gentle heat flux change [29]. By doing this, the 
characteristic of non-uniform distribution of heat flux can 
be fully considered in calculations. In this paper, the wa- 
ter wall is firstly divided into 33 sections and each sec- 
tion is further divided into equal length which is equiva- 
lent to the space step size Az. Fig. 7 shows the length of 
each section in which the heat flux is uniformly distri- 
buted. Take the space grid length about 0.1m, the loop 
can be divided into 560 segments, take two constants as 
the time step At to ensure the independence between time 
and space step. From 0 to 1 s, the time step is 0.1s. From 
1 to 500 s, the time step takes 1.0 s. 


Table 2 Initial calculating parameters 
Parameter Value 
Inlet pressure 10.16 MPa 
Inlet enthalpy 1343.1 kJ/kg 
Inlet mass flow rate 0.0653 kg/s 
Inlet pressure drop coefficient 0 
Outlet pressure drop coefficient 0 
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Fig. 8 Heat load distribution versus time 


The dynamic characteristics of furnace 


The dynamic instability of the furnace occurs when 
the heat flux density is changed for some reasons. The 
inlet and outlet of the loop will form a reciprocating 
cycle mode in which the flow rate increase or decrease 
when the dynamic instability occurs. If the resistance 
characteristics of the furnace and the fluctuation of the 
flow reach a certain resonance, then the flow fluctua- 
tions will always be in a stable state of fluctuation and 
will not decay. In some cases, the temperature of the 
tube is changed frequently, which leads to fatigue fail- 
ure of pipe and is not conducive to the safe operation 
of the boiler, so it should be avoided. In the case of 
considering heat storage of the tube metal, applying 
1.2 times heat load disturbance on the sixth loop, the 
fluctuation of the inlet flow rate and outlet flow rate 
versus time is shown in Fig.9. 

It can be seen from Fig.9, the inlet flow rate and 
outlet flow rate versus time were reversed phase fluc- 
tuation, when the inlet flow rate increases, the outlet 
flow rate decreases, when the inlet flow rate decreases, 


the outlet flow increases, the amplitude of the inlet 
flow rate and outlet flow rate decreases gradually with 
time, eventually the outlet flow rate equaled to the 
inlet flow rate and restore to a steady state, suggesting 
that the flow condition is stable. 


Analysis of the temperature 


Without the heat load disturbance, the inlet pressure of 
the sixth loop at steady state is 10.16 MPa, the inlet en- 
thalpy is 1343.097 kJ/kg and the inlet flow rate is 0.065 
kg/s. On the basis of steady state, 1.2 times heat load 
disturbance is applied on the sixth loop. The simulation 
results are shown as follows: 

Fig.10 shows the average temperature of fluid at 
steady state and the wall temperature in three cases. At 
4.1 s, the heat load reaches the maximum value and the 
inlet flow rate drops to the minimum value; but at 13.1 s, 
the inlet flow rate reaches the maximum value and the 
heat load recoveries to the initial value. In Fig.10, the 
inflection point of each curve is the phase transition point. 
Take the steady state as an example, there are two inflec- 
tion points in the curve of average temperature distribu- 
tion, one is in the 8th section and the other one is in the 
26th section, which means the fluid in the first seven 
sections is in the single phase liquid state, the fluid which 
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Fig.9 The flow versus time 
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Fig.10 The average temperature of section 


is in the section of 8~26th is in the two-phase region, the 
fluid in the remaining sections is superheated steam. 
From 1 to 500 s, the maximum temperature of the tube 
wall is 435°C, which is in the allowable range, so the 
operation of boiler is safe and reliable. 

Fig.11 and Fig.12 exhibit the temperature fluctuation 
of the 55th control volume (single phase region) and the 
275th control volume (two-phase region). From Fig.11 
and Fig.12, the fluid temperature and wall temperature 
fluctuated with time finally remained at constant values. 
The wall temperature in the single phase liquid zone is 
only 1.5°C higher than the fluid temperature, but the 
value of temperature difference between the wall temper- 
ature and the fluid in the two-phase region is 2.8°C. 
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Fig.11 The temperature of 55th segment versus time 
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Fig. 12 The temperature of 275th segment versus time 
Conclusion 


In this study, a general model considering the heat 
storage of the tube wall metal is proposed. In order to 
further expand the study of supercritical flow instability 
of boiler, a numerical program applying time-domain 
method is developed. The program is capable to simulate 
the flow instability of boiler with various geometric 
structures and thermal parameters and calculate the tem- 


perature of tube wall. Comparisons with the data of Dy- 
nastab numerical calculation show that the present model 
is suitable and adequate for the analysis of flow instabil- 
ity in water wall and the temperature calculation of the 
tube wall, so it can be used for practical engineering ap- 
plication. 

The numerical analysis of the flow instability in water 
wall of a supercritical CFB boiler with annular furnace 
show that applying 1.2 times heat load disturbance on the 
loop, the inlet flow rate and outlet flow rate versus time 
were reversed phase fluctuation and eventually the outlet 
flow rate equaled to the inlet flow rate, the wall tempera- 
ture fluctuated with time also remained at constant values, 
suggesting that the flow condition was stable. In addition, 
the results also showed that in the case of considering the 
heat storage, the flow in the water wall is easier to return 
to stable state than without considering heat storage. 
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